Manuscript submitted to Website: http://AIMscicnces.org 

AIMS' Journals 

Volume X, Number OX, XX 200X pp. X-XX 



IS THE SUBDOMINANT PART OF THE ENERGY SPECTRUM 
DUE TO DOWNSCALE ENERGY CASCADE HIDDEN IN 
QUASI-GEOSTROPHIC TURBULENCE? 



Eleftherios Gkioulekas 

Department of Mathematics 
University of Central Florida 

Box 161364 
Orlando, FL 32816-1364, USA 

Ka Kit Tung 

Department of Applied Mathematics 
University of Washington 
Box 352420 
Seattle, WA, 98195-2420, USA 



Abstract. In systems governing two-dimensional turbulence, surface quasi- 
geostrophic turbulence, (more generally cv-turbulence), two-layer quasi-geostrophic 
turbulence, etc., there often exist two conservative quadratic quantities, one 
"energy" -like and one "enstrophy"-like. In a finite inertial range there are in 
general two spectral fluxes, one associated with each conserved quantity. We 
derive here an inequality comparing the relative magnitudes of the "energy" 
and "enstrophy" fluxes for finite or infinitesimal dissipations, and for hyper or 
hypo viscosities. When this inequality is satisfied, as is the case of 2D turbu- 
lence, where the energy flux contribution to the energy spectrum is small, the 
subdominant part will be effectively hidden. In sQG turbulence, it is shown 
that the opposite is true: the downscale energy flux becomes the dominant 
contribution to the energy spectrum. A combination of these two behaviors 
appears to be the case in 2-layer QG turbulence, depending on the baroclinicity 
of the system. 

1. Introduction. The characteristic feature of two-dimensional turbulence is that 
there are two conserved quantities, kinetic energy and enstrophy. This led Kraich- 
nan [37], Leith [39], and Batchelor [4] to conjecture that there will exist two inertial 
ranges, one located upscale of the spectral region of injection and another on the 
downscale side of injection. In the upscale side, it is assumed that there is only an 
upscale flux of energy, and no flux of enstrophy. On the downscale side, likewise, 
there is only a downscale flux of enstrophy, and no flux of energy. One then uses a 
dimensional analysis argument to calculate the energy spectrum E(k), where it is 
assumed that in each inertial range E(k) depends only on the corresponding single 
flux and the wavenumber k. The same type of argument was used in the energy 
cascade of three-dimensional turbulence [3, 33, 34] . Although three-dimensional tur- 
bulence also has two conserved quantities, energy and helicity, one has the option 
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to inject energy without injecting helicity. In two-dimensional turbulence it is not 
possible to inject energy without injecting enstrophy and vice versa, because the 
two quantities are related. 

Initial efforts to simulate the enstrophy cascade yielded confusing reports of var- 
ious numerical slopes. Consequently, alternative theories have been proposed over 
the past 30 years to explain them [48, 54, 55]. Recently, in carefully set up simula- 
tions, it was shown that it is possible to obtain the enstrophy cascade in agreement 
with the KLB theory [32,44,51]. A numerical simulation with very good diagnos- 
tics has shown that the inverse energy cascade can be obtained accordingly [6]. 
There are also however many papers that question the universality of these results 
[14-16,63,64]. A review can be found in Tabeling [61] and Gkioulekas and Tung 
[27]. 

Further confusion has resulted from efforts to explain the observed energy spec- 
trum of the atmosphere with the KLB theory. Observations show that there is a 
robust energy spectrum with slope —3 which transitions at large wavenumbers into 
slope —5/3 [21, 22, 49, 50]. In the KLB theory, on the other hand, one expects that 
at small wavenumbers the energy spectrum will have slope —5/3 from the inverse 
energy cascade, which will then transition at the forcing wavenumber, into a —3 
slope from the direct enstrophy cascade. The apparent contradiction between these 
two predictions has led to various explanations and debate [17,40-43,70]. 

It was conjectured by Tung and Orlando [67] that the observed atmospheric en- 
ergy spectrum results from the downscale cascade of enstrophy and energy injected 
at the large scales by baroclinic instability and dissipated at the smallest length 
scales. If r) uv is the downscale enstrophy flux and e uv is the downscale energy flux, 
it was suggested that there is a downscale energy cascade and a downscale enstro- 
phy cascade that coexist on the downscale side of injection. Then, their separate 
contributions to the energy spectrum would give the latter a compound spectral 
shape, with a —3 slope transitioning to a shallower —5/3 slope as the wavenum- 
ber increases. The transition from —3 slope to —5/3 slope occurs at the transition 
wavenumber kt- The order of magnitude of kt can be estimated by dimensional 
analysis and it is given by k t « \/r] uv /e uv . 

General Circulation models have been shown to be capable of reproducing the 
Nastrom-Gage spectrum in agreement with observations [35,36,59]. Although the 
nature of the nonlinear interactions which give rise to the downscale energy flux 
changes from quasi-geostrophic to stratified three-dimensional in the mesoscales, as 
far as the energy spectrum E(k) is concerned it is the existence of a downscale energy 
flux which gives rise to the k~ 5 ' 3 slope, regardless of the character of the motion. 
The recent interest, typified by Lindborg [42], in understanding the fc~ 5 / 3 slope in 
terms of three-dimensional stratified turbulence is justified, since it is necessary to 
account for length scales less than 100km in wavelength where the quasi-geostrophic 
assumption fails and vortex stretching becomes important. It is the view of the 
authors that it is equally important to understand why the quasi-geostrophic model 
is capable of supporting a downscale energy cascade with A: -5 / 3 scaling, because one 
also has to account for the transition of the spectral slope from —3 to —5/3, which 
occurs between 1000km and 700km, and the upper part of the —5/3 slope below 
700km. These are presumably still within the range where quasi-geostrophy is valid. 

Tung and Orlando [67] have demonstrated numerically that a two-layer quasi- 
geostrophic channel model with thermal forcing, Ekman damping, and hyperdif- 
fusion can also reproduce the fc~ 3 — > fc -5 / 3 compound spectrum. The resolution 
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of these simulations goes down to 100km in wavelength. The diagnostic shown in 
figure 7 of Tung and Orlando [67], shows both the constant downscale energy and 
enstrophy fluxes coexisting in the same inertial range. Recent measurements and 
data analysis by Cho and Lindborg [13] have confirmed the existence of a downscale 
energy flux and estimate r/ uv ~ 2 x 10~ 15 s~ 3 and e uv w 6 x 10~ n km 2 s~ 3 . From 
these estimates we find the mean value of the transition scale k t — r] uv /e uv w 
0.57 x 10~ 2 km _1 and X t = 2ir/k t ~ 1 x 10 3 km which has the correct order of 
magnitude. 

This theory of a combined downscale enstrophy cascade and downscale energy 
cascade is contrary to the widely accepted misconception that the argument by 
Fj0rt0ft [20] forbids a downscale energy flux in two-dimensional turbulence al- 
together, and through the isomorphism theorem of Charney [12] also in quasi- 
geostrophic turbulence. Various aspects of this misconception have been clarified 
by Merilees and Warn [47], Tung and Welch [69] and Gkioulckas and Tung [24]. 

As has been pointed out by previous authors [7,18], as long as the dissipation 
terms at large-scale and small scales have finite viscosity coefficients, and the in- 
ertial ranges exist, the downscale enstrophy flux will be accompanied by a small 
downscale energy flux, and the upscale energy flux will be accompanied by a small 
upscale enstrophy flux. Dimensional analysis arguments are premised on the as- 
sumption that these additional fluxes can be ignored, consequently the energy spec- 
trum predictions obtained by such arguments are valid only to leading order. While 
ignoring sublcading effects can be justified for strictly two-dimensional turbulence, 
we will argue in this paper that for models of quasi-geostrophic turbulence, such 
as the two-layer model, the subleading contributions can be important in the in- 
ertial range and cannot be safely ignored. Predicting the form of these subleading 
corrections requires a subtle mathematical argument that goes beyond dimensional 
analysis. 

Such an argument was given previously by the authors [25, 26] as follows: the 
subleading fluxes are associated with a subleading downscale energy cascade and 
a subleading inverse enstrophy cascade that contribute linearly to the total energy 
spectrum in addition to the dominant contributions. To establish this association, 
we begin with the observation that the underlying MSR theory [46, 52] of two- 
dimensional turbulence is a linear theory and allows two homogeneous solutions 
that can be linearly superimposed: a solution that represents a downscale enstrophy 
cascade and a solution that represents a downscale energy cascade. Then, we argue 
that the existence of a small downscale energy flux, when the viscosity v is finite, 
implies the presence of a linear combination of both solutions in inertial range. 
Furthermore, due to the linearity of the MSR theory, the two homogeneous solutions 
are independent of each other, so the downscale energy cascade is independent of the 
downscale enstrophy flux r] uv and the downscale enstrophy cascade is independent 
of the downscale energy flux e uv . 

It follows that, in the downscale inertial range, the total energy spectrum E(k) 
has the following three contributions: 

E(k) = E<$ (k) + E$ (k) + E%> (k), VMo » 1, (1.1) 

with^(fc), E ( 2\k) given by 

^«»W =fl ra £ m^ 5 ^^uv(k^uv) (1-2) 

E ( 2\k) = b uvV H 3 k- 3 [ X + m(Ho)r 1/3 ©W(HW). (1.3) 
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Here Eu V '(k), E%/(k) are the contributions of the downscale energy and enstrophy 
cascade, with Duv and D^J describing the dissipative corrections. We use the 
logarithmic correction of Kraichnan [38], adjusted by the constant \ of Bowman 
[8] for the contribution of the enstrophy cascade. We have also assumed without 
explicit justification that we may ignore the possibility of intermittency corrections 
to the subleading downscale energy cascade. For the downscale enstrophy cascade 
intermittency corrections have been ruled out by Eyink [19]. For the downscale 
energy cascade we conjecture that intermittency corrections are small for the same 
reasons as in three dimensional turbulence. The scales £uv , iuv are the dissipation 
length scales for the downscale energy and enstrophy cascade. Finally, E^J(k) is 
the contribution from the effect of forcing and the sweeping interactions. The latter 
can become significant via the violation of statistical homogeneity caused by the 
boundary conditions (see Ref. [23] for details). 

Thus, in the inertial range, if it exists, where the effect of forcing and dissipation 
can be ignored, the energy spectrum will take the simple form in the downscale 
range: 



We see that the energy spectrum will take the slope of —3 for small k , and —5/3 
for large k . The transition from one slope to the other occurs at k t , given by 



In two-dimensional turbulence, the fluxes e uv and r] uv are constrained by an 
inequality that was communicated to us by Sergey Danilov [26]. This constraint 
implies that the contribution of the downscale energy cascade to the energy spec- 
trum is overwhelmed by the contribution of the downscale enstrophy cascade and 
cannot be seen visually on a plot. This result was conjectured earlier by Smith 
[60] who claimed that the downscale energy cascade can never have enough flux to 
move the transition wavenumber kt into the inertial range. The two-layer model is a 
different dynamical system than the two-dimensional Navier-Stokes equations, and 
the validity of the Danilov inequality in the two-layer model is not obvious [26, 66]. 

In the present paper we will show that in the two-layer quasi-geostrophic model, 
when the Ekman dissipation coefficient ve is below a critical value, then the Danilov 
inequality will be satisfied. We will also argue that the asymmetric presence of 
Ekman damping on the bottom layer but not the top layer may cause the violation 
of the Danilov inequality for larger values of the Ekman dissipation coefficient. In 
this case, the top layer has more enstrophy than the bottom layer, as is realistic 
in the atmosphere, and provided that the difference in enstrophy between the two 
layers is large enough, the downscale energy cascade will be made observable in the 
energy spectrum. We derive specific conditions on how large this difference needs 
to be in order for the Danilov inequality to be violated for some wavenumbers k. 
The simulation of Tung and Orlando [67] has shown that it is possible to have 
an observable downscale energy cascade, which implies a violation of the Danilov 
inequality. The role of the argument in this paper is to explain how and why this 
can happen, given that it is a surprising and very unexpected result. We also hope 
that our paper will raise renewed interest in understanding the phenomenology of 
the two-layer quasi-geostrophic model. 

In general, the exposure of the downscale energy cascade requires that the dynam- 
ics of the underlying model must be baroclinic. The role of asymmetric dissipation 
in inducing baroclinicity was discussed previously by Arbic [1,2]. We will show that 



E(k) « a^/3 fc -5/3 + buvV l/ v 3 k ~ 3 { X + ln(M„)]- 1/3 . 



(1.4) 
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in the surface quasi-geostrophic model, which represents the most extreme case of 
baroclinicity, the downscale energy cascade becomes completely dominant. An im- 
mediate implication of our argument is that the existence of an extensive observable 
m ener gy spectrum of the atmosphere has the physical interpretation that 
the atmosphere is very baroclinic. 

It should be noted that the main results of the present paper rely on the statement 
that the transition wavenumber k t from —3 scaling to —5/3 scaling is given by 
h ~ \/r]uv/£uv This statement follows from the superposition principle, but it 
doesn't require it and can also be established with the closure model of L'vov and 
Nazarenko [45]. Furthermore, the dimensional estimate k t ~ y/r) uv /e uv for the 
transition wavenumber k t has been confirmed by Tung and Orlando [67]. Possible 
criticisms of the superposition principle were discussed in a previous paper [27], 
but they are not relevant to the argument of the present paper. Extending the 
superposition principle to the two-layer QG model (the main focus of this paper), 
which has three conservation laws instead of just two, remains an open theoretical 
question. 

The paper is organized as follows. The Danilov inequality is reviewed in section 
2 where we make some simple generalizations. Its implications for two-dimensional 
turbulence are discussed in section 3. The surface quasi-geostrophic model is dis- 
cussed in section 4 and the two-layer model in section 5. Conclusions and some 
further remarks are given in section 6. 



2. The Danilov inequality in one-layer models. The governing equation for 
a wide range of one-layer hydrodynamic models takes the form: 

^ + J(V,,C)=D + J, (2.1) 

where J is the forcing and D is the dissipation and ( = —Lip. Here, £ is a linear 
isotropic operator involving the derivatives with respect to the horizontal coordi- 
nates. For a general combination of hyper- and hypo-diffusion: 

D = -^o(-A)fC-^i(-A)- /l C, (2.2) 

with p, h, positive integers, p = l,h = yields the combination of molecular 
viscosity and Ekman damping. 

For 2D turbulence, £ is given by £ = —A, where A is the Laplacian operator 
and the streamfunction tp is related to the 2D nondivergent velocity as 

For barotropic QG turbulence, also known as Charney-Hasegawa-Nima (CHM) tur- 
bulence [11, 28, 29], £ is given instead by £ = — A + A 2 , where A 2 is a given positive 
constant. Another interesting family of one-layer models arc the a-turbulcncc mod- 
els where £ = A" with A = (—A) 1 / 2 . The case a — 1 corresponds to surface 
quasi-geostrophic turbulence (SQG) which is an extreme baroclinic model, and not 
a barotropic model like 2D turbulence or CHM turbulence. 

2.1. Conservation laws. Let ((/)) be the spatial and ensemble average of f(x,y) 
defined as 

«/))= // (f(x,y)) dxdy. (2.4) 
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There are two inviscid quadratic invariants for (2.1), which are: 

A = (l/2)«(-^C)», (2-5) 
B = (1/2)((C 2 )). (2.6) 

Note that B is conserved for all linear operators L, whereas the conservation law of 
A requires that L be self-adjoint, i.e. 

((f (fig))) = «(£/)<?». (2.7) 
For example in 2D turbulence it is well-known that, 

E^(l/2)(((u 2 +v 2 ))) = (l/2)((\Vyj\ 2 )) (2.8) 

= (1/2)«HK)»=4 (2-9) 
is the kinetic energy of the 2D fluid, and 

G^(l/2)((C 2 »^i?, (2.10) 

is the enstrophy. In section 4 we see that A and B have a different physical inter- 
pretation in SQG turbulence. 

The spectra A(k) and B(k) of the conserved quadratics A and B are defined as 

m = lA«-^c< fc », (2.H) 

S(fc)-^(((C <fc ) 2 )), (2.12) 

with ip <k and £ <fc the streamfunction and vorticity fields with all the Fourier 
wavenumbers greater than k in magnitude filtered out. The relationship, ( = —Zip, 
translates into the spectral relationships in the Fourier space 

C(k) = L(|k|)$(k), B(k) = L{k)A{k), (2.13) 

where 

$(k) = J J ^(x)e- 4k - x dx. (2.14) 

We will assume that L(k) > 0, so that both A(k) and B(k) are positive. Further- 
more, we will assume that L(k) is a monotonically increasing function of k. In 
2D turbulence, L(k) — k 2 ; in CHM turbulence, L(k) = k 2 + A 2 ; in a-turbulence, 
L(k) = fc a , and in SQG, L(k) = k. 

Furthermore, A(k) and B(k) satisfy the following spectral equations: 

oMk ) + S jum_ DA(k) + FA(k) 



dt dk 
dB(k) dn B (k) 
dt + dk 



= -D B (k) + F B {k). (2.16) 



Here, DA(k) and D B (k) are the spectral dissipations rate of A(k) and B(k), respec- 
tively, with 

D B (k) = L(k)D A (k), (2.17) 
D A (k) = [v Q k 2p + Vl k- 2h ]A(k) > 0, (2.18) 

for a combination of hyper- and hypo- viscosities. Furthermore, F A (k) and F B (k) 
are the spectra of forcing also related by F B (k) = L(k)F A (k), and, H A (k) and 
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ns(fc) are the spectral fluxes of A and B. The Leith [39] constraint on the fluxes 
generalizes to 

The conservation laws for A and B arise from the following boundary conditions on 
U A (k) and IL B (k): 

lim U A (k) = lim U A (k) = (2.20) 

lim n B (fc) = lim n B (fc) = 0. (2.21) 

2.2. The Danilov inequality. Assuming that the injection (forcing) of A and B 
occurs in [ki,k 2 ], then at steady state, we have, from (2.15) and (2.16): 

r + oo 

D A (q) dq, for k > k 2 (2.22) 

Ik 

r+oo 

U B (k) = J D B (q) dq, for k > k 2 (2.23) 

U A (k) = - [ D A (q) dq, for < k < k x (2.24) 
Jo 

U B (k) = - ( D B (q) dq, for < k < k u (2.25) 



n A (k) = f 

Jk 



since F A (k) = and F B (k) — for < k < k\ and k > k 2 . For wavenumbers 
k > k 2 , we have therefore 

r + oo 

L(k)U A (k)-U B (k)= / [L(k)-L(q)]D A (q)dq<0. (2.26) 

J k 

Similarly, for wavenumbers < k < k\, we have: 

L(k)U A (k) - U B (k) = - [ [L(k) - L(q)}D A (q) dq < 0. (2.27) 

Jo 

Consequently, for all wavenumbers fee (0, ki) U (k 2 , +oo) not in the forcing range, 
we have: 

L(k)U A (k)-Il B (k) <0. (2.28) 
This inequality was brought to our attention by Danilov (2004, personal communi- 
cation) for the case of 2D turbulence. We shall see later in this paper that although 
this inequality has a trivial mathematical derivation in one-layer models, it be- 
comes nontrivial for the two-layer QG model. In all cases, what is interesting is 
the physical understanding that follows from studying this inequality, and not the 
mathematical derivations per se. 

Previously, Fj0rt0ft [20] and Eyink [18] derived a similar, but looser, bound, 
for the downscale energy flux Lf^ in two-dimensional turbulence: U.E(k) < ijo/k 2 , 
involving the total rate of enstrophy injection 770. This looser inequality is often 
used to show [58] that in two-dimensional turbulence with an infinite downscale 
range, the energy flux H B (k), vanishes. For the case of small but finite viscosity 
where the downscale spectral range is finite, the energy flux does not vanish. 

The significance of the inequality (2.28) is that it decides whether the transition 
wavenumber k t is within the inertial range, thus making a transition from the 
leading cascade to the sublcading cascade observable in the energy spectrum E(k). 
Whether this happens depends on the baroclinicity of the system, as we will show 
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below by considering three different cases: two-dimensional turbulence which is 
completely barotropic, SQG turbulence which is completely baroclinic, and the 
two-layer QG model which lies in between with respect to baroclinicity. 

2.3. Integral form of the Danilov inequality. The Danilov inequality can be 
rewritten in an integral form, by noting first that 

n fl ,„ «, . L{mA (k) - jf ' °J0>n M M 

(2.29) 

where the Leith constraint is used, ft follows that 

/ ^^-H A (q)dq = L(k)U A {k)-U B {k)<0,ykG(0,k 1 )U{k 2 ,+oo). (2.30) 
Jo Oq 

Note that the validity of the inequality is not affected by the integral going through 
the forcing range as long as the endpoint wavenumber k is outside the forcing range 
[fci, fe]. The immediate implication of this inequality is that most of the "energy" 
A has to go to small wavenumbers, provided that dL{q)/dq > 0. Unlike the flawed 
argument by Fj0rt0ft [20], this argument does not rule out the possibility of a small 
downscale energy cascade. 

Likewise, it can be shown that 



s: 



-U B (q) >0, Vfce (0,fci)U(fc 2 ,+oo). (2.31) 

Ik [HQ)] 2 dq 

A detailed discussion of the integral forms of the Danilov inequality was given in a 
previous paper [24]. 

3. The case of 2D turbulence. We begin with the case of 2D turbulence in finite 
domain with finite viscosity for the infrared and ultraviolet dissipations. This is a 
generalization of the theory by Kraichnan [37] of infinite domain with infinitesimal 
dissipation. The Danilov inequality in the inertial range is given by 

k 2 U A (k) -U B (k) < 0. (3.1) 

In the inertial range on the downscale side of injection, LLa (k) w e uv , and Tl B (k) w 
rj uv , provided that we ignore a small dissipative contribution. It follows that e ut ,fc 2 < 
r\ uv for all k in this inertial range. The energy spectrum in (1.1), valid in the inertial 
range, can be rewritten to leading order, omitting the logarithmic correction: 

E(k) « C l£ 2/3fc-5/3 + C 2)7 V3fc-3 (3 2) 

C2vli 3 k- 3 ll + ^(^r) (3.3) 




- C 2 vli 3 k~ 3 , (3.4) 

where we use e uv k 2 < T] uv and the assumption C\ < Ci- This sequence of steps is 
valid asymptotically in the limit of large separation between the forcing scale and 
the dissipation scale, for wavenumbers k in the inertial range. A similar argument 
can be made for the inertial range upscale of injection where we use 11^ (k) w — e ir 
and U B (k) s=s -r\ ir . 

It follows that on the downscale side of injection the dominant cascade is the en- 
strophy cascade with E{k) <~ k~ 3 , and on the upscale side of injection the dominant 
cascade is the inverse energy cascade with E(k) ~ fc~ 5 / 3 . By "dominant" we mean 
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that even for finite Reynolds numbers the contributions of the sublcading downscale 
energy cascade and the sublcading inverse enstrophy cascade are hidden for all the 
wavenumbers k in the inertial range. A violation of the condition C\ < C2 can 
allow, in principle, a transition to fc~ 5 / 3 scaling very near the neighborhood of the 
dissipation range. The Danilov inequality by itself cannot completely rule out such 
an effect, but so far as we know, it has never been clearly observed in numerical 
simulations of the enstrophy cascade. 

This argument supports the conjecture by Smith [60] that in 2D turbulence, on 
the downscale side of injection, we have no transition to shallower scaling E(k) ~ 
£,-5/3 jj- g th er conjecture, that the same result also holds for the two-layer QG 
model, is not true in general and will be discussed later in this paper. 

It should be noted that in the foregoing arguments it is assumed that an inertial 
range exists either upscale or downscale of injection. Unlike the case of 3D tur- 
bulence, where the downscale energy cascade is very robust, it is well known that 
in 2D turbulence there are circumstances where the leading inverse energy cascade 
[14-16, 26] or the leading downscale enstrophy cascade [63-65] may fail to appear 
as expected. Some of these issues are also relevant to the case of a-turbulence 
[62]. In general, the failure of cascades is to be attributed to the absense of a suffi- 
ciently strong large-scale dissipation sink. Since the observational evidence suggests 
that cascades exist in atmospheric turbulence, we will simply assume that without 
further discussion. 



4. The case of SQG turbulence. There has been considerable confusion over 
the physical interpretation of the surface quasi-geostrophic model. Although its 
mathematical formulation is in the form of a one-layer model, it represents a three- 
dimensional system that corresponds to the baroclinic limit of the three-dimensional 
quasi-geostrophic model. Once that is taken into account, the physical interpreta- 
tion of the spectra A(k) and B{k) and the physical implications of the Danilov 
inequality have to be revised. The viewpoint that we would like to put forth in this 
section is that whereas two-dimensional turbulence is an extreme case where the 
enstrophy cascade is completely dominant, SQG turbulence is the other extreme 
case where it is the downscale energy cascade that is completely dominant. 

4.1. 3D interpretation of SQG turbulence. As derived by Charney [12], 3D 
QG flow conserves the 3D potential vorticity £, which is advected horizontally by 
the streamfunction i[>. Here, both ip and £ are 3D fields. For constant Coriolis 
parameter /, the governing conservation law for £ takes the form: 

§ + J(^,0 = 0, (4.1) 



with £ given by 



dt 



£ = A^ + ^te?W (4.2) 



po dz \N 2 dz 



where po(z) is the ambient air density, and N 2 (z) the Brunt- V'ais'al'a frequency. 
Here we have omitted the forcing and dissipation terms. The streamfunction ip is 
also linked with the potential temperature 9 via the hydrostatic relation 

e = ^ = l d X (4.3) 

N T N dz ( ' 
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The potential temperature is also governed by a 2D-like equation 

^ + J(v ,,e ) + M! w = o, (4.4) 

where w is the vertical velocity field. 

In SQG the potential vorticity £ is assumed, a priori, to be identically zero for 
z > 0. The streamfunction ip is solved from £ = Vtp = 0. With po and -/V 2 taken to 
be constants, the horizontal Fourier transform of ip(x, y, z, t) is obtained as 

V-(k,z,t) = ^ (k,i) e - |k|W/)z , (4.5) 

using the boundedness boundary condition as z — > oo. 

Most of the dynamics in this model are occurring at the surface z = 0, where 
the boundary condition of vanishing vertical velocity w applied to the potential 
temperature (9) equation leads to: 

^ + J(V,6) = D + J, (4.6) 

where now plays the role of the conserved quantity £ in (2.1). Here we have also 
introduced thermal forcing and dissipation: D = vAQ is the thermal diffusion, and 

?=Q = a B (6 -e), (4.7) 

is the thermal heating in the commonly used form of Newtonian cooling (see Ref. 
[68]) which includes a "forcing" term a£0o and the "Ekman damping" term — a#0. 
This equation is to be solved on a 2D surface z = 0. It has the same form as the 
vorticity equation for 2D turbulence (e.g. (2.1)), except that the spectral relation- 
ship between the advected quantity and the advecting field ip is given instead 

by 

Q(k, z ,t) = l^(MKt)e-^ N ">) (4.8) 

= -|k|^o(k,i)e Hk|(W)z ] (4.9) 
= -\k\MKz,t), (4.10) 

which reduces to 6(k, t) = — |k|i/>(k, t) at z = 0. Thus SQG corresponds to L(k) — 
k. 

It has been shown by Charney [12], and more generally by Tung and Orlando 
[68], that the 3D QG energy density 



» 1 



N 2 V dz 



(4.11) 



is an invariant (i.e. independent of time), when integrated over the 3D domain. £ 
is the sum of the kinetic energy density Zk and the potential energy density £p 
which are given by 

L K = (l/2)po(« 2 + v 2 ) = (l/2)p |W>| 2 (4.12) 
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For SQG, using Parseval's identity, the energies integrated over the horizontal sur- 
face are given by 



Ep = «£p» = ^o«e 2 )>, 
Ek = «£*» = ^o«|W| 2 » 



2 Po 



(ik^(k,t)) • (-ikV*(k,t))> rfk 



/('- 



fc^(M)| : 



rfk = - Po ((e 2 )) = e p . 



(4.14) 
(4.15) 
(4.16) 
(4.17) 



It is thus seen that the kinetic energy density and the available potential energy 
density, when integrated horizontally, are equipartioned, and that 



2B = ((6 2 )) = (Ep + E K )/ Po = E/po, 
is the total energy at the lower surface. The 3D energy is, instead 



(4.18) 



E$d = 



-f 

Jo 



((£))dz = / Po dz((Q 2 
Jo 

Podz J J dk x dk. 



e| z=0 



-2|k|(iV//)z 



= ~ 2 P0 



dk^dk 



f 



f 



v N\k\ 



ee* 



z=0 



= - 2 p^ (((-V>e|, 



=o)» = Po(f/N)A, 



with A defined earlier as A = (l/2)(((—ipQ))}. Previous authors have made use of 
the similarity between the form of vorticity equation (2.1) in 2D turbulence and 
the temperature equation (4.6) in SQG turbulence to identify, by analogy, A as the 
"energy" and B as the "enstrophy" [30,53]. As pointed out in Tung and Orlando 
[68] and also here, 2B is the total energy integrated over the lower surface, and 
includes kinetic plus available potential energy. The physical interpretation for A 
was not given, but can now be seen to be the total energy integrated over the 3D 
domain. There is no potential enstrophy (£ 2 /2) per se in SQG turbulence, because 
potential vorticity £ has been taken to be zero identically. Consequently the flux of 
potential enstrophy in SQG is exactly equal to zero, thus k 2 UE 2D (k) — Ha(k) > 
and k 2 U E3D (k) - U G (k) > 0. 



4.2. The energy spectrum in SQG turbulence and a-turbulence. The ar- 
gument of the previous section can be extended to the case of SQG turbulence and 
to the even wider case of a-turbulence, which includes both the case of SQG for 
a = 1, and the case of 2D for a = 2. Here L(k) = k a and we assume a > 0. Using 
the linear superposition principle [25] and assuming the existence of inertial ranges, 
the spectrum of A(k) and B(k) are, in the downscalc inertial range: 

A(k) w C 1 {e A ) 2,i k-i+^ a + C 2 (e B ) 2/3 k-?-i a (4.19) 
B(k) = \k\ a A(k). (4.20) 
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Here, ea and eb are the constant fluxes on the downscale side of the forcing range. 
The inequality (2.28) becomes k a £A < £s for all k in the inertial ranges. Conse- 
quently, for wavenumbers k in the inertial range, the spectrum A(k) is given by 

A(k) « C 2 ( £B f^k-i-^ + g (^) 2/3 ) (4-21) 

-^(es) 2 /^-^^, (4.22) 

again in the limit of large separation between the forcing scale and the dissipa- 
tion scale and provided that C\ < C 2 . Thus, in the downscale range, there is no 
observable transition and therefore: 

A(k) « C 2 (E B ) 2 / 3 k-i-i a , (4.23) 

B{k) w C 2 (£ B ) 2/3 k-^ + i a . (4.24) 

For the SQG model, in the downscale inertial range, the energy spectrum is the 
same as that for 2B{k) in a-turbulcncc with a = 1: 

E{k) = 2B(k) w Cie^V 1 + C 2 £2D fc _5/3 - (4-25) 

Here e^d = £a is the downscale energy flux of 3D energy A and e 2 d = £s is the 
downscale energy flux of 2D energy J5. From the Danilov inequality we learn that 
the visible energy spectrum in the inertial range downscale from the injection scale 
is given by 

E{k) w C 2 £ 2 2 / D 3 k- 5/3 . (4.26) 

This k~ 5 ^ 3 energy spectrum is now predicted by our theory. The flux e 2 d is not 
the "enstrophy" flux but is the 2D flux of 2D energy 2B = E at the z = layer. 
Numerical simulations have reproduced this energy cascade and appear to indicate 
a small deviation from the k~ 5 ' 3 slope due to conjectured intermittency corrections 
[10]. 

4.3. The relevance of SQG for the troposphere and the Nastrom-Gage 
spectrum. In the troposphere, because of fast vertical mixing, the potential vor- 
ticity is almost homogenized and is very close to zero. It takes a sharp jump to a 
large positive potential vorticity in the stratosphere across the tropopause. So if 
we take the tropopause to be the "surface" and take z as positive pointing down- 
ward, we have a situation approximated by SQG for those wavenumbers which 
are large compared to the distance between the tropopause and the ground. For 
these wavenumbers, which satisy k(N/f)z ^ 1 in Eq (4.8) for all z between 
and the tropopause height of H, the ground appears to be infinitely far. Since 
the SQG solution applies in this case, and since we have predicted that the energy 
spectrum has the slope of —5/3 due to downscale flux of 2D potential energy, this 
then may explain the observed Nastrom-Gage spectrum for the wavenumbers larger 
than l/(HN/f) ~ l/1000km, since H ~ 10km and N/f is typically 100 in the 
troposphere. For smaller k, the presence of the lower ground needs to be taken 
into account, and a nonlinear Eady problem needs to be solved. Two layer models, 
on the other hand, are a good approximation for the scales which "see" both the 
tropopause and the ground, and therefore are good candidates for simulating the 
—3 part of the Nastrom-Gage spectrum. 

It is important to note that when viewed this way, the presence of the —5/3 slope 
in the Nastrom-Gage spectrum is suggesting that the troposphere is very baroclinic 
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in the vertical. This situation cannot be approximated by a two-dimensional turbu- 
lence model. When using a two-layer model to approximate this baroclinic situation, 
the two layers must be very different. In the next section we argue that this can 
be achieved with asymmetric Ekman damping. A difference between SQG and a 
two-layer model is that the downscale flux in the former is finite even as viscosity 
goes to zero from above. In the latter the finite downscale energy flux needs to be 
maintained by a finite viscosity, because of the well-known absence of anomalous 
dissipation in two-layer models. 

5. The case of two-layer QG turbulence. The results in the previous sections 
demonstrate that QG turbulence can exhibit a variety of behaviors. Barotropic 
models usually possess an energy spectrum with —3 spectral slope, as in 2D turbu- 
lence, while SQG turbulence, which is baroclinic (with its exponential decay with 
height), is expected to have a spectrum with a —5/3 slope, in the downscale inertial 
range. Two-layer QG models have both a barotropic and a baroclinic component 
[56-58], and we therefore expect a mixture of —3 and —5/3 slopes depending on the 
degree of baroclinicity. In terms of the Danilov inequality, we expect that it will 
hold for small wavenumbers and fail for large wavenumbers for sufficiently baroclinic 
systems. 

A relatively realistic two-layer model applicable to studying atmospheric turbu- 
lence in the troposphere was adopted in discussed by Tung and Orlando [67]. In this 
model forcing is due to thermal heating, which injects energy directly into the baro- 
clinic part of the total energy. The two-layer fluid sits atop of an Ekman boundary 
layer near the ground, which introduces Ekman pumping in the lower layer [31] 
but not in the upper layer. If one artificially adds an identical Ekman damping in 
the upper layer it can be shown that Danilov's inequality applies, and we leave the 
proof to the interested reader. Consequently, asymmetric Ekman damping is the 
only possible mechanism that can lead to a violation of the Danilov inequality for 
quasi-geostrophic models. 

This significance of asymmetric Ekman damping on the dynamical properties 
of quasi-geostrophic models was also discovered independently by Arbic [1,2], who 
investigated a set of similar but distinct inequalities. For the more general case of 
three-dimensional models, there are additional mechanisms, which act synergisti- 
cally with asymmetric dissipation. 

For the multiple-layer QG model, the governing equations can be rewritten in 
the form of the conservation law with £ = —Hip, if we make £ into a matrix and tp 
into a column vector. We will discuss the more general theory of multi-layers in a 
future paper. Here, we only wish to explain why and how the Danilov inequality 
can fail in the two-layer QG model. 

5.1. The formal setup. Two-layer QG models conserve potential vorticity in each 
layer in the absence of forcing and damping. In the forced-dissipative case, the 
governing equations read: 



Bottom layer: 



Top layer: 




(5.2) 



(5.1) 



where 



Ci = A?/>i - -f{ipi 



fa), C2 = AV'2 + ^(V>l-^2), 
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are the potential vorticity in each layer. k R = (2\/2f)/(hN) is the Rossby radius 
of deformation wavenumber and is taken as a given constant (h is the height). The 
dissipation terms, D,, include momentum dissipation of relative vorticity, Aipi, in 
each layer, and Ekman damping from the lower boundary layer: 

Di = ^(-Af+Vi, ©2 = K-A) p+ V2 - u E Aifa. (5.3) 
The forcing terms can be shown to satisfy 

?i = +§Q, ^ = -0Q, (5.4) 

where Q is the radiative heating term in the temperature equation. 

The two inviscid quadratic invariants are the energy E and the total layer po- 
tential enstrophics G± and G 2 given by 

£=«^Ci + V2C2» (5-5) 
G\ = ((Ci))> G 2 ^((CD). (5.6) 

The energy and enstrophy spectra are defined as 

E(k) = ^«Vi <fc Ci <fc » + J^2 <fc C 2 <fc »> (5-7) 

= ^«Ci <fc Ci <fc », (5-8) 

G 2 (fc)^^((C 2 <fe C 2 <fe », (5-9) 

and the total enstrophy spectrum G(k) is G(k) = G\(k) + G 2 (fc). We define the 
streamfunction spectra 

Ui{k) = ±((tf k tf% (5.10) 

U 2 (k)^^-((^<% (5.11) 

C(k) = ^^ k )), (5.12) 

and also ?7(/s) = U\(k) + t/ 2 (fc). It is easy to show the triangle inequality 2C(k) < 
U (k). The energy and enstrophy spectra are related with the streamfunction spectra 
via 

E(k) = {k 2 + k%/2)U(k)- k 2 R C(k) (5.13) 
G(k) = (fc 4 + k 2 k 2 R + k%/2)U{k) - k\(2k 2 - k 2 R )C(k). (5.14) 

From the standpoint of examining the Danilov inequality, it is most convenient 
to work with the streamfunction spectra. However, following Salmon [56-58], for 
physical understanding it is useful to work with the energy and enstrophy spectra. 
Furthermore, it is helpful to distinguish between barotropic energy and baroclinic 
energy as follows: Let = (ip\ + ^ 2 )/2 and r = (tpx — ^2)/2. So, tpi = ij) + r and 
?A 2 = ip — t. Now we define three spectra Ex(k), Ep(k), and Ec(k) in terms of ip 



HIDDEN DOWNSCALE ENERGY CASCADE IN QG TURBULENCE 



15 



and t: 

E K (k)^2k 2 ^U< k ^< k )) 1 (5.15) 
E P (k)^2(k 2 + k 2 R )^-((r< k r< k )), (5.16) 

E c (k)^2k 2 ±-U< k T< k )). (5.17) 

Here E^(k) is the barotropic energy spectrum and Ep(k) the baroclinic energy 
spectrum. It is easy to show that the definitions are self-consistent, i.e. E(k) = 
Ex(k) + Ep(k). The relation between the energy spectra and the enstrophy spectra 
can be now written in terms of two constraint equations: 

Gi(fc) = (l/2)[k 2 E K (k) + (k 2 + k 2 R )Ep(k)} + (k 2 + k 2 R )E c {k) (5.18) 

G 2 (k) = (l/2)[k 2 E K {k) + (k 2 + k 2 R )Ep{k)] - (k 2 + k 2 R )E c (k). (5.19) 

Thus we see that the physical interpretation of Ec{k) is that it represents the 
difference in potential enstrophy distribution between the two layers, and it is given 

by 

G 1 (k)-G 2 (k ) 
2(k 2 + k 2 R ) 

5.2. Controlled necessary condition. The right-hand-side of the Danilov in- 
equality reads 

r-+oc r+oo 

k 2 U E (k)-U G (k)= dq{k 2 D E {q)-D G {q)) = dq A(k,q), (5.21) 

Jk Jk 

so the key question is whether A(fc, q) is positive or negative for wavenumbers 
k < q < k max . Here fc max is either the truncation wavenumber in the numerical 
model, or, in the theoretical case of infinite resolutions, is the hyperviscosity dissi- 
pation wavenumber, beyond which the spectral enstrophy dissipation rate becomes 
negligible. 

With some straightforward but tedious calculations, it can be shown that the 
dissipation rate D E (k) for the total energy, and the dissipation rate Dc(k) for the 
total potential enstrophy are given by 

D E (k) = 2vk 2p+2 U(k) + 2v E k 2 U 2 (k) (5.22) 

D G (k) = 2vk 2p+2 E{k) + v E k 2 [(2k 2 + k 2 R )U 2 {k) - k 2 R C(k)}, (5.23) 

thus it follows that 

A(k,q) = k 2 D E (q)-D G (q) (5.24) 

= 2vq 2 ^ 2 [k 2 U(q) - E(q)} + v E k 2 R q 2 C(q) (5.25) 

+ v E q 2 (2k 2 -2q 2 -k 2 R )U 2 (q) (5.26) 

= 2vq 2 P+ 2 {(k 2 - q 2 )U(q) + (k 2 R /2)(2C(q) - U{q))] (5.27) 

+ v E k 2 R q 2 {C(q) - U 2 (q)) + 2v E q 2 {k 2 - q 2 )U 2 (q). (5.28) 

The first and third terms in this expression are always negative. Consequently, a 
necessary condition for violating the Danilov inequality is that the second term has 
to be positive, i.e. C(q) — U 2 (q) > 0. Otherwise, if C(q) — U 2 {q) < 0, then the 
Danilov inequality will be satisfied. A physical interpretation of this condition will 
be given in section 5.4. 
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A controlled sufficient condition to satisfy the Danilov inequality can be derived 
in terms of the physical parameters of the problem by noting that C(q) — U 2 (q) < 



(l/2)(C/i(g) - U 2 (<?))• It follows that 

A(fc, q) = 2vq 2 ^ 2 [{k 2 - q 2 )U(q) + (k 2 R /2)(2C(q) - U(q))] (5.29) 

+ u E k 2 R q 2 {C{q) - U 2 (q)) + 2v E q 2 (k 2 - q 2 )U 2 (q) (5.30) 

< 2vq 2 *+ 2 \{k 2 - q 2 )U{q) + (k 2 R /2)(2C(q) - U(q))\ (5.31) 
+ u E k 2 R q 2 {1/2)^^) - U 2 (q)) + 2v E q 2 {k 2 - q 2 )U 2 {q) (5.32) 

= \2vq 2 ^ 2 (k 2 - q 2 ) + (l/2)^4g 2 ][/ 1 (g) (5.33) 

+ [2{vq 2 P+ 2 + v E q 2 )(k 2 - q 2 ) - {l/2)u E k 2 R q 2 ]U 2 {q) (5.34) 

< \2vq 2 ^ 2 (k 2 - q 2 ) + {l/2)u E k 2 R q 2 ]UM- (5-35) 



Here, we have used the inequality 2C(k) < U (k) to eliminate the (k R /2)(2C(q) — 
U(q)) term. We have also eliminated the U 2 {q) term because it is uncondition- 
ally negative. This leads to the following controlled sufficient condition to satisfy 
Danilov's inequality: 



v E < Avk 2 ^ 




(5.36) 



Equivalently, a necessary condition to violate Danilov's inequality is 



v E > 4f fc^ax 




(5.37) 



It is interesting to note that in the numerical simulation of the two-layer model 
the algorithm adopted by Tung and Orlando [67] for determining the magnitude of 
the hyperviscosity coefficient is v E 3> ^fcnf ax , for all but the last twenty wavenumbers 
k in the dissipation range. Tung and Orlando [67] obtained an energy spectrum 
with the compound slope configuration and the transition wavenumber fc t occured 
in the incrtial range downscale from injection in agreement with the condition k t w 
\ZVuv/^uv, thus implying a violation of Danilov's inequality. 

5.3. An uncontrolled necessary and sufficient condition. The question now 
arises: is it possible to derive a sufficient condition to violate the Danilov inequality 
of the form v E k 2 R > Kvk 2 ^^ 2 for some universal constant A? So far as we know, 
this is not possible. However, it is possible to derive an uncontrolled necessary and 
sufficient condition for violating the Danilov inequality. 
We begin with defining 

2C(q) = X(q)U(q) (5.38) 
U 2 (q) = u(q)U(q) (5.39) 

Here < u(q) < 1 and —1 < X(q) < 1. We may thus rewrite everything in terms of 
U(q) by employing 



C(q)-U 2 (q) = (\(q)/2-u(q))U(q) 
2C(q) - U{q) = (\(q) - l)U(q) 



(5.40) 
(5.41) 
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Then we can rewrite A(k,q), as follows: 

A(k, q) = 2vq 2 ^ 2 \{k 2 - q 2 )U(q) + (k 2 R /2)(2C(q) - U(q))} (5.42) 

+ v E k 2 R q 2 {C{q) - U 2 (q)) + 2v E q 2 (k 2 - q 2 )U 2 {q) (5.43) 

= q 2 U(q)[-vq 2 P(2(q 2 - k 2 ) + 4(1 - \{q)) (5.44) 

+ v E k R (X(q)/2 - u(q) + 2{k/k R ) 2 u{q) - 2{q/k R ) 2 u{q))]. (5.45) 

It easy to see that a sufficient condition to get A(k,q) > is 

^4(A( q )/2- u ( g )+2(fc/M 2 «(9)-2(g/M 2 «(9))>^ 2p (2(g 2 -fc 2 )+4(i-A(g))). 

(5.46) 

The necessary condition to violate the Danilov inequality C(q) — U2(q) > 0, which 
was derived previously implies that X(q)/2 — u(q) > 0. Provided that we assume 
the stronger condition 

X(q)/2 - u(q) + 2(k/k R ) 2 u(q) - 2(q/k R ) 2 u(q) > 0, (5.47) 

we may rewrite our sufficient condition as: 

VEk 2 R 2(q 2 -k 2 ) + k 2 R (l-X(q)) 

vq 2 P ~ X{q)/2 - u(q) + 2(k/k R ) 2 u{q) - 2{q/k R ) 2 u{q) ' [ °' ° J 

To violate the Danilov inequality at wavenumber k, this sufficient condition must 
hold for all q such that k < q < k max . Since the numerator is always positive, the 
condition (5.47) is in fact a stronger necessary condition for violating the Danilov 
inequality. 



5.4. Physical interpretation of necessary conditions. We would like now to 
discuss the plausibility of the necessary condition X{q)/2 — u(q) > and the stronger 
necessary condition (5.47). To this end, we rewrite these conditions equivalently in 
terms of the physical energy and enstrophy spectra. 

It is easy to write U\(k), U2(k), and C(k) in terms of Ex{k), Ep(k), and Ec(k): 

Ui(k) = ^{{^ + r)< k ^ + r)< k )) (5.49) 

= ^ <? ^ <fe » + 2 i^ <V<fc )> + i« r<V<fe » ( 5 ' 5 °) 
_ E K {k) E P (k) E c (k) 

2fc 2 + 2(P + A:|) + k 2 ' [ ' 

and with a similar argument we find 

rr (k) E ^ | E C (k) 

U ^ k > ~ + 2WT^j (5 ' 5 ) 
c(!:) _ E K (k) E P (k) 

C[k >- 2P 2(k 2 + k%y (5 ' 53) 
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We may thus write A(fc) and u{k) in terms of Ex{k), Ep{k), and Ec(k): 

E K {k) _ E P {k) 
2C(k) 2k- 2(k 2 + k 2 R ) 

X{k >~ U{k) ~ 2 E K (k) + E P (k) (5 ' 54) 



and 



fc 2 ' (k 2 + k 2 R ) 
(k 2 + k 2 R )E K (k)-k 2 E P (k) 
{k 2 + k 2 R )E K {k) + k 2 E P {ky 



E K {k) E P (k) E c {k) 



(5.55) 
(5.56) 



(5.58) 



u(k) _U2(k) 2k 2 2(fc 2 + fc|) k 2 

u[k) -~m~ E K (k) Ep(k) - (5 - 57) 

k 2 (fc 2 + fc 2 ) 
1 (fc 2 + k 2 R )E K {k) + k 2 E P (k) - 2(fc 2 + k 2 R )E c {k) 
~ 2 (fc 2 + k 2 R )E K (k) + k 2 E P (k) ■ 

and the necessary condition A(g)/2 — u(q) > can now be rewritten as 

U wo r\ (I 2 + k R )Ec(q)-q 2 Ep(q) 

A(?)/2 - ^ = (q 2 + k 2 R )E K{ q) + q 2 Ep { q) > °" ( 5 ' 59 ) 

The denominator is obviously positive, consequently the condition is equivalent to 
Gi(«) - G 2 (q) - 2(q 2 + k 2 R )E c (q) > 2q 2 E P (q), (5.60) 

for all q such as fc < q < k max . 

This is a very interesting result. The requirement, in part, is that there should 
be more enstrophy on the top layer than the bottom layer, i.e. G\(q) > G2(q)- It is 
reasonable to expect this if there is more dissipation on the bottom layer than the 
top layer. This is the case for the model we are considering where there is Ekman 
damping at the bottom layer but not at the top layer. However, the actual condition 
also requires that the difference should be larger than the potential energy Ep(q) 
multiplied with 2q 2 . In our model, as has been pointed out by Salmon [56-58] , the 
energy is injected into the system at large scales as baroclinic energy, and most of 
this energy is converted into barotropic energy near the Rossby wavenumber kp. It 
is therefore reasonable to expect that this necessary condition will be satisfied for 
wavenumbers q > kp. 

It should be noted that this is only a necessary condition. A sufficient condition 
would require furthermore that inequality (5.48) be valid. The necessary precondi- 
tion for for the inequality to have the desired direction is the stronger requirement 



u(q) \ k R/ \ k R 

This condition can be rewritten equivalently as 

G 1 (q)-G 2 (q)-2q 2 E P (q) ^ fqY Jk 



2G 2 (q) + k 2 R {E K (q) - E P (q) ) \kp J \kp 



> 2 hr- - 2 h- • ( 5 - 62 ) 
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and it can be simplified further to give 

- (l + 4^pr^G 2 (<?) > 2q 2 E K (q) - 2k 2 (E K (q) - E P (q)). (5.63) 

The constraint u(q) > implies that 

2G 2 (q) + k 2 R (E K (q) - E P (q)) > 0, (5.64) 

and subtracting this inequality from the necessary condition above gives the follow- 
ing simplification: 

G 1 (q) - (1 + 4(q/k R f)G 2 (q) > 2q 2 E K (q), (5.65) 

for all q such that k < q < k max . 

This condition, which is also a necessary condition for violating the Danilov in- 
equality places an even stronger constraint on the difference between the enstrophy 
between the two layers, which is dependent on the ratio 4(g/fc^) 2 . In the numerical 
simulations of Tung and Orlando [67], the variability of A(q/kR) 2 is relatively small 
because they only go down to 100km in wavelength. It should be noted, of course, 
that the quasi-geostrophic model is not valid at length scales much smaller than 
100km in wavelength. When three-dimensional effects become relevant, the conser- 
vation of enstrophy is violated, and additional means for continuing the downscale 
energy cascade to smaller scales come into play. 

6. Conclusions and Discussion. The classical KLB theory of 2D turbulence 
relies for its mathematical simplicity and elegance on two unrealistic assumptions: 
that the domain is infinite, and that the Reynolds number approaches infinity. 
When these two assumptions are relaxed, the situation becomes more complicated. 
The downscale enstrophy cascade is accompanied with a hidden downscale energy 
cascade, and similarly the inverse energy cascade is accompanied with a hidden 
inverse enstrophy cascade. This is true as long as the leading cascades themselves 
exist, which requires the presence of a sufficiently strong dissipation sink at small 
wavenumbers. The fluxes associated with the subleading cascades are constrained 
by the Danilov inequality, and as a result the subleading cascades cannot contribute 
large enough terms to the energy spectrum to create an observable effect. This 
picture represents a generalization of the KLB theory to finite inertial ranges and 
finite dissipations. This situation changes, however, in baroclinic models of quasi- 
geostrophic turbulence. 

The surface quasi-geostrophic model represents an extreme baroclinic case where 
the entire behavior in the three-dimensional domain is constrained by the behavior 
of the system at the z = layer. In this model there is no enstrophy, and the 
dominant feature is the downscale energy cascade. 

We have shown that in the two-layer quasi-geostrophic model, the violation of 
the Danilov inequality is possible only as a result of asymmetric Ekman damping 
operating on only one of the two-layers. This creates an imbalance between the 
amount of enstrophy accumulated in one layer versus the amount accumulated in 
the other layer, and the downscale energy cascade will become observable on the 
condition that this imbalance is sufficiently large. We have derived in the present 
paper a sufficient condition for not violating the Danilov inequality which explains 
why the fc~ 5 / 3 spectrum has not been observed in some of the previous simulations 
of the two-layer model. We have also derived a necessary and sufficient condition for 
violating the Danilov inequality, but it is an uncontrolled condition. The numerical 
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simulation by Tung and Orlando [67] has confirmed that a double cascade with 
the transition wavenumber located in the inertial range can be realized. This can 
only occur when the Danilov inequality is violated for some wavenumbers k in the 
inertial range. The parameterization of the Ekman damping in that simulation does 
in fact satisfy the necessary condition derived in this paper. 

As long as we operate within the framework of multiple- layer models with a finite 
number of layers, one cannot rule out the alternative theory that the atmospheric 
energy spectrum might reflect a double downscale cascade of helicity and energy 
instead of enstrophy and energy (see discussion in section 6.5 of Ref. [9], and figure 
3 of [5]). However, most of the current debate has been focused on the somewhat 
mysterious nature of the very extensive and robust fc~ 5 / 3 spectrum. 

Our work in the present paper explains why it can be reproduced in numerical 
simulations that use baroclinic models, while the same effect cannot be realized in 
simulations of two-dimensional turbulence. On the other hand our work here does 
not rule out the possibility that the shallower part of the spectrum observed by 
Nastrom and Gage [49] over the mesoscales can be due to dynamics other than QG, 
whether it is barotropic or baroclinic, especially for scales of 100 km or less (see 
e.g. Ref. [42,43] with Bousinesq dynamics). Our present work serves to point out 
that over the larger scales Qt 600km), where the transition to a shallower spectrum 
occurs, baroclinic QG theory by itself is a viable mechanism for explaining the 
transition from —3 to —5/3 slopes. 

Furthermore, as proposed first by Tung and Orlando [67], the downscale energy 
flux, which is important in explaining the fc -5 / 3 energy spectrum over the mesoscales 
in most theories, originates at larger scales (the synoptic scales). Its contribution 
to the energy spectrum is hidden for smaller wavenumbers under the fc~ 3 part of 
the spectrum, and then emerges for larger k past the transition scale. It remains an 
open question, one that is beyond the scope of this paper, to understand in detail 
how this downscale energy flux can be continued into length scales too small for 
QG theory to describe, and how it is eventually dissipated. 
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